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Abstract. 

C*") ' The relativistic theory of structure formation in cosmology is based mainly on linear 

perturbations about a homogeneous background. But we are now driven to understand 
the theory of higher-order perturbations in full detail, both from observational and 
' theoretical points of view. An important aspect of this lies in defining gauge-invariant 

, perturbations at any order. We present a new covariant approach to this based on 

(— | | a local separation of tensors into their scalar, vector and tensor parts. Such a local 

decomposition necessarily requires a trio of mutually annihilating differential operators 
q , which form the basis for defining gauge-invariant objects. It makes no use of non-local 

S-^ ' Green's functions or pre-defined gauges, and can be used to define families of scalar, 

52 ■ vector and tensor modes at any order one chooses. 

c3 

> 

pr) ; 1. Introduction 

^ , A challenging problem in relativistic cosmology is perturbation theory around a 

qq . homogeneous and isotropic background when extended beyond first-order. A particular 
difficulty lies in preserving general covariance under perturbations, which necessitates 
considering general coordinate transformations order-by-order alongside a pertubative 
expansion, and trying to eliminate gauge degrees of freedom. For a tensor to be gauge- 
invariant (GI) at a given order, it must vanish (or be constant) at all lower orders (see [1] 
for a full discussion). It is widely acknowledged that using tensors which are fully gauge 
invariant is the best way to construct higher-order perturbation theory. One important 
result of Nakamura's [2] has been to show how starting from gauge-invariant first-order 
quantities in the metric approach one can generate GI objects at second- and higher- 
order. Thus one can define a second-order GI tensor perturbation of the metric in the 
Poisson gauge, for example, but this perturbation is gauge- dependent in the sense that 
it depends on using the Poisson gauge to begin with; starting in a different gauge one 
has a different tensor perturbation, even though it too can be GI (see [3] for details). 

An important restriction for these gauge-based GI quantities is that they are 
defined non-locally, which requires assumptions about unknowable boundary conditions 
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at infinity and inherently relies on assumptions about conditions outside our horizon. 
While this is fine in a mathematical sense, it is at the very least peculiar in cosmology 
because such conditions really are, in a fundamental sense, unknowable, and cannot 
influence conditions within our horizon. 

In addition to this, it is often not clear what the GI quantities mean in the following 
sense: consider a second-order GI extension to the usual Newtonian potential; because 
this is GI it must correspond to the leading term of some tensor which vanishes at first- 
order (since we know GI variables are tensors which vanish at lower orders), but we do 
not know what this object is nor how to find it. 

All constructions of higher-order perturbation theory so far are based on the so- 
called metric approach to perturbation theory, whereby the metric is split into various 
parts whose dynamics is governed by different components of the field equations. An 
important alternative approach to linear perturbation theory has been the 1+3 covariant 
approach, which has gauge-invariance built in [4,5]. But this approach has languished 
at first-order, unable to compete with the metric approach because of the difficulty of 
defining GI variables at higher-order, except in special circumstances. This problem we 
rectify here by providing a new general technique for defining GI objects at any order 
covariantly - and these quantities can be used in the metric approach too. 

This method offers the following advantages: it is conceptually simple, being based 
on permutations of three mutually annihilating differential operators; it is local as there 
are no Green's functions of \7 2 + nK involved; it is covariant, so no gauges or coordinates 
are necessary to define these quantities; and finally, variables are defined such that they 
are automatically scalar, vector or tensor modes. 

2. The covariant approach to cosmological perturbations 

In cosmology there exists physically defined observers and therefore reference frames 
with which to consider physical quantities. Any observable quantity necessarily relies 
on such an observer being chosen. This underlies the 1+3 covariant approach to GR 
initiated by Hawking, Ehlers and Ellis [4,5]. Any observer can find, say, the unique 
CMB reference frame, or the energy frame (in which the net heat flux vanishes), and 
use such a frame to define the electric and magnetic parts of the Weyl tensor which do 
not depend on any coordinate system, thereby preserving general covariance. That is, 
there are various velocity fields around which physically and invariantly define the rest 
of the spacetime quantities; the coordinate system one uses remains irrelevant. Such 
considerations have led to the 1+3 covariant approach to cosmological perturbation 
theory in which gauge-invariance and frame-invariance are subtly separated. 

The 1+3 approach is a semi-tetrad conversion of the field equations of GR into a 
set of evolution and constraint equations which are derived from the Ricci and Bianchi 
identities. The variables involved are all covariant objects (i.e., they are tensors), defined 
through projections with a physical velocity field u a and its associated spatial projection 
tensor h a b = g a b + u a Ub and volume element e a bc = u d Vabcd [4]. All tensors may then 
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be irreducibly decomposed into scalars, vectors and projected symmetric and trace- 
free (PSTF) tensors. These rank-0 tensors, projected rank-1 tensors and PSTF rank-2 
tensors collectively describe the dynamics of the spacetime and are governed by a system 
of first-order PDEs. In cosmology, where there is approximate spatial homogeneity and 
isotropy, such a projection leads to the background spacetime being described by a 
family of invariant scalars, such as the energy density and expansion rate, provided u a 
is the velocity of the fluid. 

More usefully than this is what happens at first-order. Any of the 1+3 objects with 
an index must vanish in the background as otherwise it would break the symmetry, and 
so it must be GI by the Stewart- Walker Lemma. This important theorem states that a 
tensor is GI provided it vanishes in the background [6]. Except in specialised situations 
little work has been carried out at higher order [7, 8] because of the difficulty in finding 
GI quantities. 



2.1. Local scalar-vector-tensor decomposition 

All rank-1 and -2 tensors used here are orthogonal to u a , and rank-2 tensors are 
projected with h a b, symmetric and trace-free which we denote using angle brackets 
on indices. We define the spatial covariant derivative acting on scalars or spatial tensors 
as D a X b ^ c = h a a 'h b b ' ■ ■ ■ h/V a 'X b / ... c / . The irreducible parts of the spatial derivative of 
PSTF tensors are the divergence, curl, and distortion, defined as [9] 

divA b ... c =B a X ab ... c (1) 
curl X ab ,„ c = e de{a D d X b c) e (2) 
disX ca ,„ b = D {c X a _ b) . (3) 

Then, the spatial derivative of a rank-n PSTF tensor Xa„ = X aia2 ___ an may be 
decomposed as (for n = 1, 2, 3) 

2t?/ — 1 n 

D b X An = - — — div X(An-i h an )b — curl X^a^e )b c + dis X bAn . (4) 

in + 1 n + 1 1 

Note that the divergence decreases the rank of the tensor by one, the curl preserves it, 
while the distortion increases it by one (and all are PSTF). Keeping this in mind one 
can drop the indices on differential operators as long as it's explicit the valance of the 
PSTF tensor which is being acted on. 

At maximal perturbative order a GI object may be considered as a field on an 
FLRW background. Then, a general rank-2 GI PSTF tensor can be split into a non- 
local scalar, vector and tensor parts 

X ab = S ab + V ab + T ab = D <a D 6) S + D {a V b) + T ab , (5) 

where the scalar part S ab is curl- free, the vector part V ab is solenoidal, D a V a = =>- 
D a D b V ab = D a D 6 D( a Vb) = 0, while the tensor part is transverse, divT a = 0. A similar 
decomposition exists for rank-1 tensors, but we shall concern ourselves here with rank- 
2 tensors only: a PSTF rank-2 tensor can be formed from a rank-1 one by taking a 
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distortion without affecting its SVT classification. An important problem using such 
variables is that they are non-local [10]. 

It was demonstrated in [8] that there exists equivalent local variables corresponding 
to S a b, Vab and T ab which can be found by acting with suitable combinations of differential 
operators on X ab . These are local because they do not rely on the solution to an 
elliptic differential equation that the usual SVT decomposition (5) implicitly relies on; 
consequently there is no reliance on integrals over all space and unknowable boundary 
conditions. These differential operators, which act on rank-2 PSTF tensors, are:i" 



S = dis D div div (6) 

V = dis curl div (7) 

T = [dis div - |D 2 + f curl 2 ] curl . (8) 

These are constructed so as to preserve the rank of X ab .§ Note that div T = 0, 



div div V = and curl S = as required. Then, SX ab depends only on S ab , VX ab 
depends only on V ab , TX ab depends only on T ab ; similarly, TX ab is a divergence-free 
rank-2 tensor, hence a bone fide tensor mode. These operators are mutually annihilating 
when acting on a quantity at maximal perturbative order (i.e., acting on the background) 
in the sense that 

SV = VS = ST =TS = VT = TV = 0. (12) 
3. Gauge-invariant objects 

At first-order finding gauge-invariant objects is trivial: any object with an index must 
be GI by the Stewart- Walker lemma. How do we find GI objects at second- and higher- 
order? The trick in going from zeroth to first lies in taking gradients of scalars: scalars 
are non-zero in the background but their gradients will not be, thereby providing useful 
GI objects (in addition to things like the electric and magnetic Weyl tensors which are 
GI anyway). At first-order, however, all of the standard 1+3 covariant quantities are 
non-zero, so it has been a long standing problem as to how GI second-order objects may 
be formed. Can we take derivatives of first-order quantities in a suitable way such that 
they are zero until we get to second-order? What about higher-orders? 

One method is to excite only a certain degree of freedom at first-order, such as 
scalar modes [7,8]. Then, any variable which is a pure vector or tensor mode has to be 
second-order and, hence, gauge- invariant. For example, in this case we can easily see 

| Note that T looks a bit different from in [8] . We have substituted for the curvature in terms of curl 2 
and D 2 so that T can be used unambiguously as a differential operator in any spacetime. 
§ Written out in full these are: 

SX ab = B {a B b) B c B d X cd (9) 
VX ab = e cd(a B b) B c B e X e d (10) 
TX ab = \e cde B (a B e B c X b) d + \e cd{a [ZB b) B e B c X d - B 2 B c X b) d - AB c B 2 X b) d 

- B e B c B e X b) d + B e B c B b) X e d + B c B e B b) X e d + B e B b) B c X e d + B c B e B d X b)e ] (11) 
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that the vorticity is a pure vector as it is divergence free - consequently it must be GI 
at second-order; similarly, VE ab is also a GI vector mode at second-order, and TH ab is 
a GI tensor mode. What is not obvious is how to isolate the scalar modes which are 
induced at second-order in a GI way. This actually illustrates the main problem: how 
do we handle the general situation when all modes are excited at first? 

The key lies in analysing what happens in going from zeroth to first. Consider the 
gradient of the energy density D a p, which is GI at first-order. Despite appearances, 
D a p is actually a mixture of a scalar and a vector mode. We can see it contains a 
vector mode because if we take the curl, we find curl D a p oc u a . Similarly, to isolate the 
scalar part of D a p we have to take a divergence. The reason for this is that because p 
itself is not GI when we take its gradient, the covariant derivative is not a background 
covariant derivative but has first-order connection terms mixing in; the vector degree of 
freedom arises from this. To summarise, we must first take a derivative orthogonal to 
the symmetry of the background to form a GI object; then we must differentiate again 
in two ways to isolate the scalar and vector part. 

Now let us generalise this to one order higher. Consider disoj a ;,: this is a pure 
vector degree of freedom at first-order because u a is divergence free. Consequently, at 
first-order, we must have Sdisu ab = and Tdisu ab = 0. Therefore, at second-order 
SdisUab and Tdisu ab must be GI PSTF tensors. Thinking of S, V and T as 'orthogonal 
operators', taking derivatives orthogonal to V (since dis u a b is a pure vector) results in 
GI variables at the next perturbative order. Given either of these objects we can operate 
once more with S, V or T to isolate pure scalar, vector or tensor degrees of freedom. 

More generally, if we take any PSTF rank-2 tensor such as the shear or electric 
Weyl curvature, we can operate with S, V and T to form a set of 3 new PSTF objects: 
in the case of the shear for example we have So~ ab , Vo- a b an d To~ a b- At first-order these 
have an invariant meaning in terms of the SVT decomposition, but at second they do 
not. Now operate on each of these three with S, V and T, to form a total of 9 new 
rank-2 tensors. Of these, SSa ab , Wcr ab and TTa ab are not much use to us as they are 
non-zero at first-order. However, all the mixed cases have to vanish at first-order and 
therefore must be GI at second. That is, an operator to produce a GI PSTF tensor at 
second-order from any PSTF first-order one belongs to the annihilating set 

n = {sv,sr,vs,vr,rs,TV} (13) 

After this operation, they still do not have an invariant meaning in terms of the SVT 
split because S, V and T must act on GI quantities for this to be the case (analogously 
to D a p not being a scalar at first-order). To extract the SVT parts of an element of 
n[cr a6 ], say, we must differentiate again with S, V or T, giving a total of 6 scalars, 6 
vectors and 6 tensors. Explicitly, we have: 

scalar: SIl A a ab -w SSV(r ab , SSTa ab , SVSa ab , SVTa ab , STSa ab , STVa ab 
vector: VIl A a ab *** VSVo- ab ,VSTa a b,VVSa ab ,VVTa a b,VTSa ab ,VTVo- ab 
tensor: TIl A a ab -w TSVa ab ,TSTa ab ,TVSa ab ,TVTa ab ,TTSa ab ,TTVa ab , 
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where A — 1, ... ,6, forming a set of 18 GI objects which now have an invariant meaning 
in terms of the SVT split. Similar quantities can be formed starting from any PSTF 
rank-2 tensor, including things like D( a D^p. 

Let us now extend this to third-order. Again we can form objects which vanish 
at second-order by applying annihilating projections of S, V and T to those SVT GI 
variables just defined - i.e., SVHa, STUa, VSTIa, VTIPa, TSTIa, TVHa give a total 
of 6 x 6 = 36 objects which vanish at second-order and consequently must be GI at 
third. Note that we have essentially acted with II twice on a first-order PSTF tensor in 
order to find all these objects. Further action with S, V or T will then split these into 
sets of scalars, vectors and tensors. Thus, 7~VSVSo- a b is a GI third-order tensor mode. 

We can generalise this to any order. We have seen that at second-order II acting 
on a first-order PSTF tensor gives a set of 6 GI variables, while at third we have to act 
with II§n, i.e., the outer product of II with itself. In index notation we have ILaIIb. 
At perturbative order n we must act with II n — 1 times, or 

n ®(n-i) = n <g> n <g> • • • <g> n . (14) 

S V ' 

n-i times 

In index notation we have instead IF^IL^ • • ■ IIa„_ 1 , giving a total of 6 n_1 possible GI 
PSTF tensors. The SVT parts can be separated with a final action of S, V and T as 
before. 

4. Conclusions 

We have demonstrated for the first time how to generate covariant, locally defined, 
GI objects at any order when perturbing away from an FLRW background. Utilising 
local differential operators which mutually annihilate on the background we can form 
rank-2 tensors which represent only scalar, vector or tensor degrees of freedom, provided 
they act upon a GI object at a given perturbative order. Recursively acting with these 
operators in the way described above then provides a neat way to define a large variety 
of GI variables at any given order. 

How would such GI quantities be used in practise? In principle, for a given order 
n one can calculate evolution equations for II <g) ^ n ~ 1 - ) X a fe, where X a b stands for a set of 
appropriate rank-2 tensors, using standard commutation relations and the 1+3 covariant 
equations. Products of quantities which appear will be of perturbative order n—1 and 
can be written in terms of products of H^ n ~ 2 ' X a f>; these will have their own evolution 
equations. The final system will yield a large system of PDEs which can be solved order 
by order. This is left for future work. 

The physical information these variables contain will be very subtle to interpret of 
course. When considering something like TSVo- a b f° r example, we cannot read this as 
the 'the tensor part of the scalar part of the vector part of the shear at second-order', 
because once one moves up a perturbative order the meaning of the SVT decomposition 
disappears. Analogously, Va a b is the vector part of the shear at first-order, but not 
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at second. But such difficulties are to be expected: there are only a limited number 
of tensors that are regularly used in practise which we have intuition for, so to expect 
one to 'appear' at third- or tenth-order representing a GI tensor mode which also has a 
simple physical interpretation is implausible. Nevertheless, as higher-order perturbation 
theory is developed in full generality this will be an important avenue to explore. 
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